Application of the Feshbach-resonance management 
to a tightly confined Bose-Einstein condensate 
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We study suppression of the collapse and stabilization of matter-wave solitons by means of time- 
periodic modulation of the effective nonlinearity, using the nonpolynomial Schrodinger equation 
(NPSE) for BEC trapped in a tight cigar-shaped potential. By means of systematic simulations, 
a stability region is identified in the plane of the modulation amplitude and frequency. In the 
low-frequency regime, solitons feature chaotic evolution, although they remain robust objects. 



I. INTRODUCTION 

Dilute atomic Bose-Einstein condensates (BECs) are 
accurately described by the Gross-Pitaevskii equation 
(GPE), alias the cubic nonlinear Schrodinger equation 
(NLSE) Q. The sign of the cubic term in the GPE cor- 
responds to the self-defocusing or focusing, if interactions 
between atoms in the condensate are characterized, re- 
spectively, by the positive or negative s-wave scattering 
length. The self-focusing GPE in any dimension (ID, 2D, 
or 3D) gives rise to soliton solutions, which are stable in 
the ID case. The creation of ID matter-wave solitons has 
been reported in experimental works 0, while 2D and 
3D solitons are unstable against the critical and super- 
critical collapse, respectively (these 2D states are usually 
called Townes solitons, TSs) It was predicted that 
TSs may be stabilized in the framework of the 2D GPE, 
without using an external potential, if the constant scat- 
tering length is replaced by a time-dependent one, that 
periodically changes its sign Q. In BEC, this can be im- 
plemented by means of the Feshbach-resonance manage- 
ment (FRM), i.e., by applying a low- frequency ac mag- 
netic field which acts via the Feshbach resonance [5] . This 
stabilization mechanism was demonstrated in optics, in 
terms of the transmission of a light beam through a bulk 
medium composed of layers with alternating signs of the 
Kerr nonlinearity Q , and then in the framework of the 
2D GPE [7], [lfj ■ A somewhat similar technique was pro- 
posed recently, making use of a linear coupling, induced 
by means of a resonant electromagnetic wave, between 
two different hyperfine states of atoms, which feature op- 
posite signs of the scattering length [ll| . The analysis of 
the FRM was extended to include averaging techniques 
[IH , generation of solitons from periodic waves [l3[ , the 
stabilization of higher-order solitons [l4[ , management of 
discrete arrays [151 ] , and the case of a chirped modulation 
frequency However, the stabilization based on the 
FRM may be, strictly speaking, a transient dynamical 
regime, as extremely long simulations suggest that the 
FRM-stabilized TS may be subject to a very slow decay 
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The stabilization of 3D solitons by means of the FRM 
technique alone is not possible, but stable 3D solitons 
were predicted in a model combining the FRM and a 
ID periodic potential 0]. Similarly, the stabilization is 
possible when the FRM is applied in combination with 
a parabolic potential which strongly confines the con- 
densate in one direction Most relevant to the ex- 
periment is the "cigar-shaped" setting, with the BEC 
tightly confined in two transverse directions, while the 
third direction remains free Q. In the usual approxi- 
mation, with the cubic nonlinearity in the corresponding 
ID GPE, the analysis of the FRM in the latter setting 
amounts to that reported in Refs. However, if the 
density of the condensate is not very low, the descrip- 
tion in terms of the cubic nonlinearity is inappropriate, 
the respective ID equation taking the form of the non- 
polynomial Schrodinger equation (NPSE). In particular, 
it admits the onset of the collapse in the self-attractive 
condensate in the framework of the ID description [l7l ]. 
Accordingly, a relevant problem, which is the subject of 
the present work, is to study the possibility of the col- 
lapse suppression by means of the FRM technique in the 
framework of the ID NPSE. It is relevant to mention that 
the NPSE was recently used to describe Faraday waves 
generated in the cigar-shaped trap by a time-periodic 
modulation of the strength of the transverse confinement 
[TH - We introduce the model in Section 2, and report 
results obtained by means of systematic numerical simu- 
lations in Section 3. 



II. THE NONPOLYNOMIAL SCHRODINGER 
EQUATION 

The normalized form of the 3D GPE with the trans- 
verse harmonic trapping potential, which acts in the 
{x, y) plane, is 

id t ^ = [-(l/2)(d 2 x +dl + d 2 z ) 

+ {l/2)pl + V{z) + 2vg\vb\ 2 }^. (1) 
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Here ip(x,y, z,t) is the mean-field wave function, and 
p\ = x 2 + y 2 . Further, g = 2a s N/a±is the nonlin- 
earity strength, with a s the s-wave scattering length, 
N the total number of atoms in the condensate, and 
aj_ = y/h/(mu}_i) the confinement radius imposed by 
the transverse harmonic potential of frequency lu± , with 
m the atomic mass. In Eq. ([1]) length and time are 
measured in units of a± and uij 1 - As usual, g > and 
g < correspond to the repulsion and attraction between 
atoms in the BEC, respectively, and V(z) is a weak ax- 
ial potential, which may be present in addition to the 
strong transverse confinement. Being interested in the 
stabilization mechanism that does not require the extra 
potential, we set V(z) = 0. Then, the 3D e quat ion can 
be reduced to the NPSE by means of ansatz 
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where ID wave function f(z, t) is subject to the nor- 
malization condition, f^°° \f(z, t)\ 2 dz — 1. Following 
Refs. one can eliminate the transverse width, 

a 4 = 1 + g\f\ 2 , arriving at the NPSE, 
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In the case of g = const, stationary solutions are looked 
for as / = exp (— ipyt) $(2), where /i is the chemical po- 
tential, and real function <&(z) obeys equation 
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Some numerical methods for simulations of the GPE and 
NPSE (with g = const) were presented in Ref. [22| . 

In the case of the attractive nonlinearity, g < 0, the 
form of Eq. implies that the amplitude of the wave 
function is limited from above by a critical value, 



\f\ 2 < 



f\ 



(5) 



A dynamical collapse sets in, with transverse width a 
shrinking to zero and the solution developing a singular- 
ity in finite time, as |/| 2 approaches the critical value p^ |. 



In Ref. 21(, this was called two-dimensional primary col- 



lapse, as it is related to the transverse 2D dynamics. 

In addition to the dynamical collapse, the NPSE also 
admits a static collapse, in the framework of stationary 
equation (@|: for g < and V{z) = 0, this equation ad- 
mits bright-soliton solutions only below the critical value 
of the nonlinearity strength, \g\ < \g c \ = 4/3 At 
\g\ < I (7 C I , the axial density |$(z)| 2 in the bright-soliton 
solution is smaller than the critical value imposed by con- 
dition (O. In Ref. [2l|, this kind of the collapse was 
called three-dimensional primary collapse, as it involves 
a quasi-spherical 3D soliton. With regard to the defini- 
tion of g, this restriction determines the largest number of 



atoms possible in the soliton, N < N max — (2a±/3 \a s \). 
With a± and \a s \ on the order of /xm and 0.1 nm, re- 
spectively, which is typical for experiments in the 7 Li 
condensate @, one has N max ~ 10 4 atoms. 

The FRM technique which makes it possible to sta- 
bilize 2D matter-wave solitons is based on the respec- 
tive GPE, id t i> = [~(l/2)(d 2 x + d 2 )+27rg(tM 2 ]^, 
where constant g < is replaced by a periodic func- 
tion, g(t) = g dc + sin {wt) , with \g Ac \ < g &c , so that 
g(t) alternates between attraction and repulsion. The 
stabilization requires the presence of the constant ( "dc" ) 
component which corresponds to the self-attraction on 
the average, i.e., gdc < 0. 

The action of the FRM within the framework of the 
NPSE was not considered before. To explore this situ- 
ation, we take g(t) as indicated above, arriving at the 
following modification of Eq. Q : 
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Our objective is to identify a region in parameter space 
{Saci Sdo ^} where the solitons subjected to the "man- 
agement" represent stable solutions of Eq. ©. 



III. RESULTS 

Localized solutions to Eq. ^ were categorized as sta- 
ble ones if, in direct simulations, they featured persistent 
pulsations, avoiding collapse or decay up to t = 5000 (in 
some cases, the stability was checked up to t = 80000). 
However, the application of this criterion to the case of 
u! < 0.5 is complicated by the fact that, under the action 
of the low-frequency management, the soliton tends to 
develop an apparently chaotic behavior, although with- 
out a trend to decay, see below. Fixing the time interval 
as that comprising a large number of periods, simulations 
become increasingly more difficult for lu — * 0. 

The simulations we performed by means of the Crank- 
Nicolson algorithm with open-ended boundary condi- 
tions. The initial state was fo(z) = sech(z), that was 
refined by the integration of Eq. ^ with g ac = in 
imaginary time [22j . The so generated configuration was 
then used as the input to simulate Eq. ^ in real time. 

Typical examples of stable and unstable solutions are 
shown in Figs. [T]and[2] (in the model with g ac = 0, the 
respective soliton is stable). In the latter gradu- 
ally growing amplitude of the soliton attains critical value 
([5]) at finite t, which implies the onset of the dynamical 
collapse. It happens when the argument of the square 
root in Eq. ([B]) becomes zero, i.e., the transverse width 
of the soliton shrinks to zero. Note that, with g = g(t), 
critical density ([5]) of the axial wave function does not 
necessarily correspond to the maximum of g(t). Indeed, 
in Fig. [DJa) the collapse happens at g(t) smaller than its 
maximum value, |<7dc| + |<7ac|- 
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FIG. 1: A typical example of a stable soliton solution to 
Eq. (JBJ. (a) The evolution of the soliton's amplitude in time, 
(b) A snapshot of the soliton at t = 100. The integration 
step is At = 0.05, and the size of the integration domain is 
L = 200. Parameters are <?dc = —1, Sac = 1-5, u = 1. 
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FIG. 3: (Color online). Stability borders in the plane of the 
time-modulation parameters, (o),(; ac ), as obtained from sys- 
tematic simulations for different fixed values of g^c- In cases 
when the threshold depends upon the integration time (see 
Fig. [4}, the respective symbol corresponds to the mean value, 
with the error bars given as per the respective semi-dispersion. 



FIG. 2: A typical example of the collapsing soliton. (a) The 
evolution of the amplitude up to f ~ 25, when it reaches 
critical value (J5j). (b) A snapshot of the soliton just before 
the onset of the collapse. Parameters are the same as in Fig. 
[TJ except for g ac = 2. 



Results of systematic simulations are summarized in 
the form of the stability diagram displayed in Fig. [3] 
The stability thresholds shown in the figure, i.e. the max- 
imum value of (7 ac admitting stable solitons, were found 
by slowly increasing <? ac in steps of Ag ac = 0.1, until 
the instability was attained. The shape of the stabil- 
ity domain in the plane of the management parameters, 
(<7 ac ,w), is roughly similar to that which was found in 
management models of a different type, with the time- 
periodic modulation applied not to the nonlinearity, but 
to the optical-lattice potential, which is necessary for the 
existence of stable solitons in those cases. These include 
the ID model for gap solitons, with a positive scattering 
length 23], and the 2D GPE with a negative scattering 
length and ID or 2D periodic potential, that stabilizes 
TSs in the respective settings [24]]. As in those works, 
one may expect that here, at very large values of lu, the 
stability region will start to expand in the direction of 
larger values of g ac , as in the limit of lu — ► oo the ac term 
averages to zero. 

The stability borders in Fig. [3] are not extended to 
lu < 0.5, as in the region of the low-frequency modula- 
tion the solitons feature persistent but apparently chaotic 
evolution. In fact, the stability domain is well defined 
for oj > 1, while in the intermediate region, 0.5 < to < 1, 
the randomness of the soliton evolution makes the stabil- 
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FIG. 4: Dependence of g ac at the soliton's instability border 
on the integration time, for different modulation frequencies 
and (jdc = —1- Symbols denote the stability limits at dif- 
ferent frequencies: uj — 0.5 (diamonds), = 0.75 (crosses), 1 
(triangles), 1.5 (squares). 



ity border dependent on the integration time - see Fig. 
31 which demonstrates a natural trend to a decrease of 
the effective instability threshold with the increase of the 
evolution time, if the threshold is sensitive to it at all. 



IV. CONCLUSION 

We have used the NPSE, i.e., the ID mean-field equa- 
tion for tightly trapped BEC, with the nonpolynomial 
nonlinearity admitting the onset of the collapse in the 
framework of the ID description, for the study of the 
stabilization of solitons by means of the FRM (Feshbach- 
resonant-management) technique. The results were re- 
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ported in the form of stability diagrams in the plane 
of the management parameters, (g ac ,Lu). The stability 
domain is roughly similar to that reported in linear- 
management models for ID gap solitons and 2D TSs 
(Townes solitons), supported by optical lattices subjected 



to the time-periodic modulation. However, stability do- 
mains of such a form have not been reported before in 
models of the nonlinearity management. At small values 
of the modulation frequency, the stability border becomes 
fuzzy, as solitons feature chaotic evolution in that case. 
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